Positron production scenarios and the angular profile of the galactic center 511-keV line 
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The observed angular profile of the 511-keV photon excess from the Milky Way galactic center 
can allow us to select among combinations of various dark matter and other positron production 
mechanisms with various models for the dark matter distribution. We find that a relic decay scenario 
gives too flat an angular distribution for any dark matter distribution in our survey, but that a dark 
matter-dark matter coUisional scenario, or a scenario that involves particles emitted from a localized 
central source producing positrons some distance out, can match the observed galactic center angular 
profile if the dark matter distribution is neither too flat nor too cuspy. Additionally, positron migration 
or diffusion before annihilation broadens the angular profile to an extent that an average migration of 
more than half a kiloparsec is not viable with most dark matter distributions. The observed angular 
profile is also consistent with the occurrence of transient events in the past, followed by isotropic 
positron diffusion. 



I. INTRODUCTION AND PRELIMINARY RESULTS 

The center of our galaxy is observed to be a concen- 
trated source of 511-keV gamma-rays. The source re- 
gion appears coincident with the galactic bulge, with a 
near circular symmetry with about a 6° FWHM or of 
scale about 1 kpc in FWHM diameter and with a to- 
tal flux from the bulge of about 0.89 x 10"'' photons 
cm~^ s"^ (ll0]- The identity of the positron source re- 
sponsible for these gamma-rays remains undetermined. 
There have been a number of astrophysical suggestions; 
some of these suggestions are included in In ad- 

dition, the possibility that dark matter may create these 
positrons has been widely discussed ll8l-ll2[|. 

In this paper, we focus on the dark matter possibilities 
and consider what the observed angular profile of the 
511-keV gamma-rays tell us about possible dark matter 
production mechanisms, and further consider whether 
the dark matter production scenario can inform us re- 
garding the dark matter distribution near the center of 
the galaxy. 

For purposes of this paper we think of the surviv- 
ing dark matter particles as WIMPs and will categorize 
dark matter production as being either from decay of 
relic long lived excited WIMPs j^ lTsll , from collisions of 
WIMPs where the cross section time velocity has no ve- 
locity dependence, or from collisions of WIMPs where 
the cross section has noticeable velocity dependence. 
The relic decay idea is that there are excited WIMPs 
with a small mass gap from the ground state WIMP and 
whichproduce positrons when they decay to the lower 
state l9]- The lifetime of the excited WIMP needs to be 
long enough to allow many of them to survive since 
freeze-out, yet short enough there would be enough de- 
cays at present to account for the observed rate, even 
after considering that only a fraction of the surviving 
WIMPs would be the excited ones. 

The collisional mechanism depends on having final 



states that lead directly or by some chain to positrons. 
One example is WIMP pairs armihilating to lighter dark 
matter gauge boson pairs, where the dark bosons have 
a weak coupling to the visible sector so that they soon 
enough decay into e+e" pairs, XX ~^ 'P'P ~^ e~^e~e~^e~ . 
This in only an example. For velocity independent col- 
lisional case the spatial distribution of positron produc- 
tion does not depend on the reaction specifics. We only 
need that the cross section and dark matter density suf- 
fice for obtaining the overall rate. 

We also consider the possibility that there is a lo- 
calized source of initiating particles that spread out 
and at some distance interact in reactions that pro- 
duce positrons. Photons streaming outward clearly 
bremsstrahlung positron-electron pairs, but one can 
imagine dark matter alternatives. As a possibility, if 
WIMPs connect to the luminous world via gauge par- 
ticles ("dark photons"), one could produce e+e" pairs 
via a Bethe-Heitler-like process where the WIMP plays 
the role of the nucleon and the dark gauge boson plays 
the role of the virtual photon that couples to the nu- 
cleon fT3]. The example is just an example; the angu- 
lar profile generated by any mechanism where a central 
source produces particles that then produce positrons 
some distance out will be the same. Since the visible 
galactic bulge is somewhat oblate, seeing a circularly 
symmetric positron hot spot in the center of the galaxy 
can be a clue that the dark matter halo is taking a strong 
role in the production mechanism. 

That the angular (or spatial) profile of the 511-keV 
gamma-rays can discriminate among dark matter mech- 
anisms is known lITsl, [Tsl - fT^ . The simple reason is that 
the total number of decays per unit volume in the de- 
cay mechanism depends on the dark matter density lin- 
early, while the number of collisions per unit volume 
depends on the same density quadratically so that as 
the observational line of sight gets closer to the center of 
a dark matter concentration, the collisional mechanism 
will peak more sharply that the decay mechanism. 
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In formulas, for the relic decay scenario the angular 
flux is given by the line of sight integral 



f*T 
'in 



dl n{r) , 



(1) 



where flux <J> is the number of 511-keV photons per cm^ 
per sec, dO/dO is per unit solid angle of source, n{r) is 
the number density of WIMPs, T is the decay rate of the 
excited WIMPs, and /* is the fraction of WIMPs in the 
excited state. (The formula is for the case of each decay 
yielding one positron.) Similarly for collisions. 



do 

dn 



In 



dl n{r) 



(2) 



for the case that av is constant and ng+ is the number of 
e+ produced per collision. 

In the scenario with a central source, the isotropic 
flux of initiating particles falls like the inverse distance- 
squared from the source. We will at least pro tern place 
the source at the galactic center, so the Irne-of-sight in- 
tegral is like the relic decay case with n{r) replaced by 
n{r)/r^. This too will give and angular profile steeper 
than the relic decay scenario. 

Additionally, before plotting any curves the angular 
flux results are smeared with a 3° FWHM gaussian since 
the data we compare to comes from a telescope with a 
3° spatial resolution jHH. 

The difference between the angular profiles from 
positrons produced by relic decays and by velocity in- 
dependent collisions is seen in the heavy dashed (colli- 
sion) and dotted (decay) curves in Fig.[ll using the NFW 
dark matter distribution model We will consider in 
the next section a selection of other dark matter distri- 
butions and see how the results differ among them. 

Also on Fig. [T] is a heavy solid line for a velocity de- 
pendent collision possibility which we shall return to 
shortly, and two light dotted lines representing the data 
from the spectrometer SPI on the INTEGRAL (Interna- 
tional Gamma Ray Astrophysics Laboratory) satellite. 
The central bulge photons have a approximate circular 
symmetry and there is also a broad 511-keV photon sig- 
nal from the galactic disk. The photon distribution from 
the bulge can be represented by a single gaussian of 
FWHM 6°, or somewhat better, by two concentric gaus- 
sians of FWHM 2.1° and 8.0°, as |2] 



give a idea of the level of the (albeit not circularly sym- 
metric) disk contribution. The normalization of the data 
curves is absolute while the normalization of the calcu- 
lated curves is arbitrary; we are interested here in the 
angular shape. 

One sees already that the dark relic positron produc- 
tion curve is too flat to match observation, and this 
persists even when we consider other models for the 
dark matter distribution. Collisional positron produc- 
tion gives an angular profile that falls more steeply and 
looks closer to the data. 

The orange dot-dash curve on Fig. [T] is for a central 
source emitting particles isotropically that in turn pro- 
duce positrons at radius r. One has 



do 

dn 



dN [ n{r)_ 



87t2 dt 



(4) 



where dN /dt is the emission rate of the initiating parti- 
cles and ctq is the cross section of the reaction that pro- 
duces the positrons. 

Finally, the heavy solid curve on Fig. [T] is from a col- 
lisional process with a significant velocity dependence 
in the cross section. Typically one models the veloc- 
ity distribution of dark matter (or other) particles as a 
Maxwellian with an upper cutoff at the escape veloc- 
ity UM- As the escape velocity gets higher closer to 
the center of the galaxy, there are more high velocity 
particles found there and the rate of collisions near the 
center is enhanced for this reason as well as for reasons 
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dn 
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cm s sr 



(3) 



for e in radians, fli = (2.l7r/180) /2.355 and fli = 
(8.07t/180)/2.355. (The supplement to d suggests 
somewhat wider gaussians, but 12] states the flux in each 
gaussian.) This is shown as the lower light dotted line 
in Fig. [l] The upper light dotted line is the same but 
with a constant 0.004 photons cm 



2 s ^ sr ^ added to 



FIG. 1: Comparing positron production mechanisms and com- 
paring to data. The green dotted line is for the relic decay sce- 
nario, the red dashed line is for the velocity independent colli- 
sional scenario, the blue solid line is for the collisional scenario 
with a velocity threshold (the WIMP mass = 500 GeV af- 
fects the profile for only this case), and the orange dot-dash 
curve is for a central source emitting particles that in turn pro- 
duce positrons at radius r. The calculated curves in this plot 
use the NFW fl^ dark matter distribution. The light dotted 
curves indicate INTEGRAL data, as explained in the text. 
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of overall density. Hence the velocity dependent colli- 
sional process rises more steeply toward the center of 
the galaxy that the velocity independent one, and the 
match to the data is better. 

The specific velocity dependent model we use fol- 
lows |@I where the velocity dependence occurs because 
of a threshold. Velocity dependence can also occur be- 
cause higher angular momenta are involved, and would 
lead to qualitatively similar results. The threshold oc- 
curs because, in the model, positron production oc- 
curs when an excited WIMP is produced by a WIMP- 
WIMP collision and the excited WIMP emits an electron- 
positron pair as it decays back to the groimd state. If the 
mass of the WIMP is m^^^ and the mass gap is S, then the 
threshold velocity is 



Vthresh= \l^6/m^. (5) 

The collisional cross section is modeled as 

^'"re] = C^O \J - ^l.resh ^(^rel " V thresh) , (6) 

where for WIMP velocities v and v' , Vyf,i = \v — v'\. The 
number of scatterings per time per density is given by 

l,CTv) (r) = I d^v d^v' f{v, f)f{v', r)a{v,,i)v,,i , (7) 

where f{v,f) is the WIMP velocity distribution at radius 
r. The calculated flux is modified from the previous re- 
sult to 

g = i^/d/.V)MW. (8) 

The velocity distribution is taken to be 

fiv,r) = Nir) exp{-v^/2vl,,) 9{vesc - v) (9) 

with 

^rms — 200 km/s and normalization J d^v fi^v^T^ — 
1. The normalization parameter N(r) depends on the 
distance r from the center of the galaxy because the es- 
cape velocity does. 

The escape velocity follows from noting the circular 
rotation velocity Vc of the galaxy is approximately con- 
stant out to some radius ri„„x- For the present paper we 
ignore the variation in the rotation curve near the center 
of the galaxy. Assuming there is no mass beyond r„iax, 
one obtains the escape velocity 

Vesc = ^J2[l-ln{r/rmax)] ■ (10) 

for r < r„iax- We use Vc = 220 km/ sec. The curve in 
Fig. [T] is drawn for r„,ax = 25 kpc. Recent analysis of 
stellar velocity data suggests that the constant circular 
speed holds out to 80 kpc [21]. Using Tmax = 80 kpc 
broadens the angular profile somewhat, but not dramat- 
ically 
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FIG. 2: Dark matter density models vs. distance from the galac- 
tic center. The blue dashed line is the ihdp profile, the solid red 
line is the Merritt et al. profile, the dotted black line is the NFW 
profile, the long-dash orange curve is the Cembranos and Stri- 
gari profile, and the dot-dash green line is the TF dark matter 
profile. The curves are normalized to give the same dark mat- 
ter density at our sun's location. 



We have thus far only used one dark matter distribu- 
tion model, and have tacitly (until now) assumed that 
the positrons annihilate close to where they are pro- 
duced. We discuss other dark matter distributions in 
the next section, and the effects of positron diffusion or 
migration are discussed in Sec. [nil with just the com- 
ment here that positron migration can broaden the an- 
gular distributions but not narrow them. 

Further, in Sec.lIVl we discuss possibilities that follow 
upon the galactic center being more active in the past 
than it is today. 



II. DARK MATTER DISTRIBUTION MODELS 

The results in shown in SecUused only one model for 
the dark matter density in our galaxy. Other dark mat- 
ter distributions have been suggested, and it is natural 
to enquire how using these would change the 511-keV 
angular profile already presented. 

We consider several mass density profiles, in partic- 
ular, the Merritt et flZ._profile ||22||, the Navarro-Frenk- 
White (NFW) profil e llgll , the isothermal halo dark mat- 
ter profile (ihdp) ||231,|24|] , a somewhat more singular pro- 
file used by Cembranos and Strigari |2^. here scaled 
to give the local dark matter density, and a truncated 
flat (TF) distribution from i26ll with parameters obtained 
in [27]. The latter halo distribution gives a flat rotation 
curve all the way to the inner regions of the galaxy, as is 
observed for the Milky Way at least at radii larger than 
a few tenths kiloparsec. 
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Explicitly, the poM profiles are 



Pihdp 

Pnfw 

PMerritt 
PCembranos 

Ptf 



Plocal 
Plocal 



R2 + r2 

R0(l + R0/r,)2 




l + (Re/ro)^ 
l + (r/ro)8 

^^0 1 



3/16 



r2 + fl2 



(11) 



where Rc = 2.8 kpc, = 20 kpc, a = 0.2, r_2 = 25 kpc, 
ro = 10 kpc, and a = 105 kpc. We took Rq = 8.5 kpc 
ar\d used pi^cal = (0-43 ± 0.15) GeV/cm^ from iH, with 
their xmcertainties added in quadrature. (The TP fit pa- 
rameters from translate to p{R^^) = 0.53 GeV/cm^, 
with error limits compatible with piocal given here; we 
are in any case interested in this work in the shape of 
the distribution.) The dark matter profiles give different 
amounts of dark matter near the galactic center, as can 
be seen in Fig. |2l 



A. Relic decay positron angular profiles 

The angular distributions of 511-keV photons that 
follow from the selected dark matter profiles and the 
relic WIMP decay scenario are shown in Pig. Ettop). 
The normalizations of the curves are controlled by 
the lifetime of the excited WIMP and by fraction of 
WIMPs in the excited state. By choice, our normaliza- 
tions give results that cross at 5°, and assuming each 
WIMP decay gives one positron, the lifetime/excited 
WIMP fractions for the displayed results are t / f* = 
(4.0,5.1,4.9,6.7,16.1) x 10^^ years for the ihdp, Merritt, 
NFW, Cembranos, and TF distributions, respectively. 
The light dotted curves in Fig, p] again represent data 
from the INTEGRAL satellite ElIlH- 

Four of the distributions are too flat to match the data. 
The result from the fifth distribution is possible at least 
in the inner region. However, this distribution itself 
ascribes the flatness of the rotation curve to halo mat- 
ter. The Milky Way rotation curve is flatfish in to even 
within 1/4 kpc, which is within 1.7° from our view- 
point, but usually one thinks that much of the gravita- 
tional pull and shape of the rotation curve in the inner 
region is coming from ordinary matter. In that case, car- 
rying a 1 /r2 dark matter profile down to the lowest r 
overestimates the amount of dark matter there. 
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FIG. 3: Angular distribution of 511-keV photons due to (top) 
decays of relic excited WIMPs near the galactic center, (second) 
due to positrons created through WIMP annihilation with (av) 
independent of velocity, (third) due to WIMPs being collision- 
ally excited and subsequently decaying back to the ground 
state by emitting e+e^, and (bottom) due to a point source 
of initiating particles that spread and produce positrons when 
they interact away from the source. The blue dashed line uses 
the ihdp profile, the solid red line uses the Merritt et al. profile, 
the heavy dotted uses the NFW profile, the long-dash orange 
curve uses the Cembranos and Strigari profile, and the dot- 
dash green line uses the truncated flat distribution. The fainter 
black dotted lines are a two-Gaussian representation of the 
data described in the text, one with a 0.004 photons/ (cm^sec 
sr) background included as in [90. 
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B. Velocity independent collision processes 

One can also imagine a scenario of positron produc- 
tion through WIMPs annihilation. The annihilation oc- 
curs without velocity threshold and in the simplest case 
with ((Jc) independent of velocity, hence, {av) can pull 
out of the line of sight integral. In this case, the angular 
profile for the 511 keV line photon is flatter than in the 
excitation case. 

In Fig. |3] (second panel), we show the angular dis- 
tribution of the 511 keV line for the five different 
dark matter profiles. In the plot we use {av) = 
271(1.7,0.3,0.2,0.01,0.5 x lO"*^) x lO^^O cmVs for the 
ihdp, Merritt, NFW, Cembranos, and TP distributions, 
respectively. The normalizations are chosen so that they 
cross at 2.5°. 



C. Velocity dependent collision processes 

We calculate the pairs intensity using five different 
dark matter profiles mentioned above. The WIMP 
mass is set to 77?^ = 500 GeV with the excited state 
is at i5 = 1 MeV above the ground state. Again, 
to approximate spatial response of SPI we smoothed 
the line of sight integral by 3° PWHM beam. The 
normalizations are chosen such that they cross at 3°. 
The scattering cross section of the displayed results are 
(70 = 271(3.11,0.51,0.41,0.032,0.0007) x lO^^e cm^ for 
the ihdp, Merritt, NFW, Cembranos, and TP distribu- 
tions, respectively. The results are shown in Fig. |3] (third 
panel). 



D. Steady central source, remote production scenario 

Another scenario envisions a localized source isotrop- 
ically emitting particles which produce positrons in in- 
teractions that take place some distance from the source. 
It is natural but not required that the source be coinci- 
dent with the galactic center, defined by the visible cen- 
ter of the baryonic matter or by the center of the dark 
matter halo. The outward streaming particles may al- 
low both baryonic or dark matter to participate in parti- 
cle creation, or may emphasize one over the other. The 
baryonic or dark matter possibilities in this scenario will 
give the same angular profile, to the extent that the bary- 
onic and dark matter distributions are equally cuspy 
and the source, baryonic matter distribution, and dark 
matter distribution are all concentric. 

The bottom panel of Pig. |3] shows angular profiles 
for this central source, remote production scenario in 
the concentric case, using the dark matter distributions 
given earlier. The angular profiles are quite steeply 
falling, since there is a f alloff both for the flux fac- 
tor and for the falloff of the target matter distribution. 
In this scenario, only the ihdp dark matter profile is too 



flat to match the data. Positron production off baryonic 
matter in this scenario should have an angular profile 
similar to the TP model, since this by itself gives a flat 
rotation curve, as is observed down to fraction of a kilo- 
parsec. 

We should expect that the central source has a finite if 
small extent and that the positrons will migrate at least 
a little before annihilating. Accordingly, we have fur- 
ther smoothed the profiles by a small amoimt (of scale 
about 10^^ kpc). This smoothing controls a mild singu- 
larity affecting some of the curves, but if the scale of the 
smoothing is small it does not visibly affect the curves 
shown in Pig. |3l The question of how longer positron 
migrations affect the angular profiles is discussed in the 
next major section. 

A possibility within the central source, remote pro- 
duction scenario is that the source is separated from 
the baryonic matter and/ or dark matter centers. (There 
are examples of "dual active galactic nuclei" 129] and 
of galaxies with bright off -nuclear X-ray sources js^]-) 
This leads to an asymmetric or ellipsoidal distribution of 
positrons, as has been reported from observation in |1^, 
as there will be two positron concentrations, one at the 
source and one at the matter concentration. The hot 
spots may not be of equal strength, ond one or both 
could be muted by the 3° PWHM detector resolution 
smearing. 

We show in Pig. |4] a dark matter only example of re- 
sults possible in this case. In this example, the source 
is taken to be the same distance from us as the relevant 
mass centers, but is offset 6° in the sky from the dark 
matter center. The upper panel shows the results un- 
smeared. The coordinates are centered on the mass cen- 
ter, and there is a bright spot there with a brighter spot 
at the location of the source. The lower panel shows the 
results with the 3° PWHM telescope resolution averag- 
ing. The dark matter center hot spot is largely smeared 
away. One might hypothesize within this scenario that 
the faint hot spot seen in the Ref . [19] data could be the 
location of a dark matter distribution center or of a dark 
matter clump. 

E. WIMP mass effects and extended halo effects 

The mass of the WIMP does not affect the angular pro- 
file from the dark matter relic or velocity-independent 
collisional scenario. In these cases the only effect of the 
mass is to change the number density for a given mass 
density, and this is merely a constant overall factor. 

With a velocity dependence, mass changes can affect 
the shape of the angular distribution. Por the thresh- 
old model we have considered, the mass gap & may be 
chosen to be just large enough to produce positrons, so 
that we think of as a constant. Then for lighter WIMP 
masses, the threshold velocity increases. Since at a given 
radius there is a given upper limit (the escape velocity) 
to the initial state velocities, there are radii where the 
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FIG. 4: 551 keV gamma-ray distributions created by particles 
emitted from a localized source and interacting in the interstel- 
lar medium to produce positrons, in the case that the source 
and the mass center of the medium are separated by 6° in the 
sky. The coordinates are centered on the mass distribution, 
the source is to the left. The upper panel is unsmeared, the 
lower panel includes a 3° FWHM gaussian smearing to mimic 
the resolution of the telescope of the INTEGRAL spectrometer. 
The contours are at constant steps in the log of the intensity. 



reaction cannot proceed at all. As the threshold veloc- 
ity increases, the region where the reaction can proceed 
shrinks to a small ball about the galactic center. Since 
we smear with a 3° FWHM gaussian to match the tele- 
scope resolution, there is a limit to how narrow an an- 
gular profile we can obtain, and w 200 GeV is light 
enough to give a result close to this limit for the & we 
have chosen. This result is shown by the solid line in 
Fig- El Four of the other lines in this fugure are reference 
lines which are the same as Fig. [T] The other new line in 
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FIG. 5: Comparing positron production mechanisms and data, 
but with m Y = 200 GeV. The dashed and dotted curves are the 
same as Fig. [l] and calculated curves use the NFW dark mat- 
ter distribution. The solid curve is for the excited dark matter 
coUisional mechanism and still uses Ymax = 25 kpc, while the 
dot-dash curve has rmax = 80 kpc. 



Fig. E] is the medium weight dot-dash line, which is the 
= 200 GeV result for rmax = 80 kpc. 



III. POSITRON MIGRATION 

We have so far assumed that positrons annihilate near 
where they are made. Positrons can migrate, for dis- 
tances depending on their initial kinetic energy and the 
topology of the magnetic fields. Most positrons with a 
few MeV initial kinetic energy traveling through matter 
densities typical of the galactic center will survive a few 
times 10^ years before stopping and annihilating from 
rest ll3Tll32ll . The energy loss, at these energies, is mainly 
due to ionization and collisional losses. The positrons 
freeze into the magnetic field lines, traveling in tight spi- 
rals around them and turning with the field lines as long 
as the turning radii are wider than the cyclotron radii 
(as is usually the case). The actual distance the positron 
travels from the source hence also depends on the topol- 
ogy of the magnetic field. If the field has a significant 
random component with a cell size of a parsec, as seems 
to be the case in the solar neighborhood fSS'], the few 
MeV positrons will travel about a kiloparsec before an- 
nihilating. For the galactic center, there is progress in 
imderstanding the magnetic field |33 - I37ll , but it is not 
yet a clearly understood subject. 

The positron migration, following random trajectories 
from the point of their creation spiraling along random 
field lines before finally stopping, will be treated as a 
random walk and approximated by Gaussian distribu- 
tion. A simple way to take this into account in the cal- 
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culation of gamma ray angular profile is by performing 
an additional 3-dimensional Gaussian smoothing. The 
effect of this is a flatter gamma ray angular distribution. 

In Fig. [6l we show the gamma ray angular profile 
in the collisionally excited dark matter scenario with 
Merritt dark matter profile for three different diffusion 
lengths: (0.1,0.3,0.5) kpc FWHM. The angular distribu- 
tion in the case of 0.1 kpc diffusion length is practically 
indistinguishable from the unsmoothed case in Fig. |3] 
and in the case of 0.3 kpc diffusion length we still have 
a moderately good fit. (For simplicity, the width of the 
Gaussian smoothing used in the calculation is constant, 
neglecting the possibility that the diffusion length may 
depend on the densities of matter in the surroundings, 
the positron's initial kinetic energy, and the topology of 
the magnetic field.) 
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FIG. 6: Angular distribution for Meritt dark matter profile 
in the collisionally excited dark matter scenario. Red-dashed 
line, blue-long-dashed line and orange-dotted-dashed line are 
the angular distribution corresponding to positron diffusion 
length of 0.1, 0.3, 0.5 kpc FWHM, respectively. 



In conclusion, positron migration may affect signif- 
icantly the angular profile of the 511 keV gamma ray 
from the galactic center. A kiloparsec diffusion length 
appears to give a too flat profile. Indeed, with the Mer- 
ritt dark matter distribution, a diffusion length of a half 
kiloparsec already seems to spread the angular profile 
imacceptably A large diffusion length also tends to blur 
the distinction among different dark matter profiles for 
the shape of of the gamma-ray angular distribution. 



IV. TRANSIENT GALACTIC ACTIVITY SCENARIOS 

The galactic center is currently not active and the mas- 
sive black hole that lies there radiates at about 8 or- 
ders of magnitude below the Eddington luminosity Jssll . 
There is evidence that the galactic center has been more 
active in the past; a few recent references include ll39l - 
lllll . Sporadic high luminosity events from the galactic 



center appear to have occurred in the past few hundred 
years i39l1 and there is speculation about sizable spo- 
radic activity millions of years in the past also. There 
have been considerations of how transient gala ctic cen- 
ter activity could create the positron excess 14211 . 

We will also consider how transient activity can con- 
nect to the galactic center positron excess. Suppose that 
there have been several catastrophic events in the galac- 
tic center. The catastrophic events may be either point- 
like positron injection in the center or illumination by 
neutral particles from a localized source, with a 1 / flux 
falloff, later producing positrons via interactions with 
the ambient matter. 

In the localized positron injection case, and we sup- 
pose positron injection occurs only at definite times, 
and the positrons diffuse. Diffusion proceeds as a ran- 
dom walk, and the positrons form a gaussian distribu- 
tion with a width that depends on the number of cells 
and the step size. The width, the step size (cell size of 
the random magnetic field), and the number of steps 
(time from the positron injection) form a trio where if 
we know two we can work out the third. For exam- 
ple, for the wider of the two gaussians in the data, the 
FWHM is about a kpc, and if the cell size is 1 parsec, the 
half-width is equivalent to 500 steps in a straight line 
or that number squared in the random walk. This puts 
the time since the injection at a quarter million parsecs 
converted to time, or about 800,000 years. The narrower 
gaussian is about a third as wide, so its initial injection 
would have occurred about 80, 000 years in the past. 

The spatial profile in this scenario is described by the 
sum 



do 
do 



(12) 



where d is the distance the positrons are from their 
starting point, a is the step size of the random walk, 
(7,- = y/ctia, and f/ is the time elapsed since injec- 
tion. If "F]NHM" is given in distance units, then 
the time since injection can be calculated from f, = 
F]NHM^ / [Sea In 2). Thus a localized injection can pre- 
cisely reproduce the observed angular profile of the non- 
disk part of the positron annihilation flux. 

The coefficients A, determine the total annihilation 
flux. They are determined by the number of positrons 
thermalized since injection. One naturally expects a 
larger fraction of in-flight annihilations originating from 
the younger, or narrower gaussian, positron population. 

In the second version, a transient localized source illu- 
minates the entire volume, possibly with energetic pho- 
tons, that in turn produce positrons from ambient parti- 
cles. The difference between this version and the "cen- 
tral source, remote production" scenario in Sec. lIIDl lies 
in the variability of the source. Here the galactic cen- 
ter region is illuminated over a time scale of lOOO's of 
years, i.e., over a very short time scale compared to dif- 
fusion times. Afterwards diffusion can occur, resulting 
in broadening of the positron distribution. 
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FIG. 7: Positron annihilation profile caused by a transient cen- 
tral source. The black solid curve is for centrally localized 
positron injection, over two different injection times. The red 
heavier dot-dashed curve represents positron production from 
ambient Merritt et al. dark matter triggered by a central radi- 
ation source, also at two separate injection times indicated by 
the lighter dot-dashed curves. Isotropic diffusion is applied to 
both cases. The green narrower dot-dashed curve corresponds 
to a more recent injection event, while the broader blue dot- 
dashed curve curve corresponds to an older event. The dot- 
ted curve is a representation of the data with a constant piece 
added to indicate the level of the disc background. 



Results of the two versions of this transient source sce- 
nario are shown in Fig. [T] Localized positrons injected 
at two distinct times give gaussian spatial profiles and 
hence, as shown by the solid curve in Fig. [7] match any 
gaussian representation of the data directly, with gaus- 
sian parameters determined in terms of positron diffu- 
sion characteristics. (The dotted curve is again a repre- 
sentation of the data with a constant piece added to in- 
dicate the level of the disc background.) The transient 
localized neutral particle source with remote positron 
production is shown as the heavier dot-dashed curve 
and is made from two contributions. One, shown as the 
light dot-dashed curve that is higher near the origin, has 
just the f allof f given by the flux dependence and 
the density falloff of the Merritt et al. matter distribu- 
tion, and the 3° FWHM detector resolution smearing. 
The other, shown by the light dot-dashed curve that is 
higher at wide angles, is the same except that it has a 
broader spread implemented with a 8° FWHM gaussian 
smearing. These contributions are interpreted as due to 
one very recent outburst and one outburst farther in the 
past, about 400,000 to 1.6 million years ago using a cell 
size of 1 /2 to 2 parsecs. 



V. SUMMARY AND DISCUSSION 

We have examined the angular distribution of the 511- 
keV gamma-rays from the galactic center, mainly in the 



context of dark matter initiated positron production and 
for a variety of models for the distribution of dark mat- 
ter in the inner regions of the galaxy. 

One of the putative mechanisms, decay of dark relic 
particles into final states containing positrons, is not vi- 
able in light of the flatness of its angular profile com- 
pared to the data. The degree of flatness differs for dif- 
ferent dark matter distribution model we have exam- 
ined, but in all cases the data is steeper still. 

Positrons produced from collisions rather than de- 
cay show a steeper falloff with angle, mainly because 
of the quadratic dependence on dark matter density, 
and can match the observed angular falloff in the galac- 
tic bulge. Additionally, a generic model where initi- 
ating particles steadily radiating from a central source 
produce positrons some distance out also gives a suffi- 
ciently steeply falling angular distribution. Within the 
latter scenario, if the radiation source is displaced from 
the center of the mass distribution, then it leads to an 
azimuthally asymmetric annihilation flux. 

At large angles, the dark matter scenarios generally 
fall well below the data. However, there are also defi- 
nitely existing astrophysical sources of positrons in the 
form of, for example, low mass X-ray binaries, pulsars, 
and /5+ emitting nuclei from supernova explosions 0]. 
The numbers of supernova produced (6+ emitters, in 
particular, are about the right magnitude to match the 
large angle (and not circularly symmetric) data. 

We also considered transient activity within the galac- 
tic center, characterized by several catastrophic events 
in the past follow by relatively quiet periods. Using 
a random walk model for positron propagation we re- 
late the spatial distribution parameters of 511-keV emis- 
sions to characteristics of the random components of the 
magnetic fields in the galactic center environment. Un- 
der this scenario, the observed gaussian spread allows 
us to determine products of the random field cell size 
and the time elapsed since the injection events. We fur- 
ther considered the possibility that a transient localized 
source triggered remote production of positrons which 
then diffuse. The resulting 511-keV annihilation profile 
is broadened by an amount dependent on the time since 
the injection event. 

Within the context of the steady state scenarios, we 
considered relic decay, collisional, and localized source 
scenarios where the positrons annihilated close to where 
the were formed and also considered effects of isotropic 
positron diffusion upon the angular profile. (We did 
not consider long range or anisotropic positron trans- 
port, which others have suggested as mechanisms to al- 
ter the correlation between production and annihilation 
regions 0-111, 113] •) A kiloparsec diffusion length ap- 
pears to give too flat a profile. For the Merritt et al. ll22ll 
or the NFW lITsIl dark matter profiles, even a half kilo- 
parsec diffusion length appears to lead to angular dis- 
tributions that could not be reconciled with observation. 
In addition, a large diffusion length tends to blur the dis- 
tinction among different dark matter profiles regarding 
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